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Abstract 

We explore approximations to the effective 5 factor for the 7 Be(p, 7) 8 B re- 
action using different approximation for the integral over the Gamov peak. 
In the temperature range of interest for solar neutrino production, S c s(E) 
may be determined to within 0.5% from 5(20) with no direct information 
on the derivatives of S(E). 
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The 7 Be(p, 7) 8 B reaction, at energies of approximately 20 keV, plays an important 
role in the production of solar neutrinos |]J . The subsequent decay of the 8 B is the source 
of the high energy neutrinos to which many solar neutrino detectors are sensitive. The 
cross section for this reaction is conventionally expressed in terms of the S factor which 
is defined in terms of the cross section, a, by 

S(E) = a{E)E exp [2irr](E)} , (1) 



where r](E) = Z\Z20L^j [ic 2 /2E is the Sommerfeld parameter, a is the fine structure 
constant, and /i is the reduced mass. The reaction rate per particle pair is || 
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where T is the temperature, k is the Boltzmann factor and b = (2n) 1 ^ 2 ire 2 ZiZ2/h comes 
from the Coulomb penetration factor. In order to study the approximations to the 
integral in Eq. @ it is convenient to introduce an effective S factor defined by || 
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where r = 3E /kT and E = (bkT/2) 2 / 3 is the location at which the function 
peaks. Since the integrand is strongly peaked the integral may be 
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done by the saddle point method 
S eS _s(E) = S(E ){l 



]. To first order in 1/r the result is [| 
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The approximation depends on S and its first two derivatives, all evaluated at E = E . 
By expanding S about E — MeV we obtain an expression where we only need S and 
its derivatives at the origin. This results in the approximation || 
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Eqs. (|]) and (|5J) both explicitly depend on the first and second derivatives of the S factor. 
Another approximation may be made in which the derivatives in Eq. ([|) are neglected, 
giving a much simpler form: 



S'eff-Snd = S(E ) < 1 + 



12r 



Finally we may replace S(Eq) with 5(20) in this last equation to yield 
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Se ff -s2o = 5(20) 1 + 



12r 
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To evaluate these approximations we need a functional form for S(E). For E < 
100 keV, the functional form of S for the 7 Be(p, 7) reaction is well approximated by [[| 

0409 

S(E)/S(0) = Q " 375 l E + 0.703 + 0.343E (8) 

where E is in MeV. With this expression for S(E)/S(0) and taking S(0) = 19.0 eVb 
H, we have calculated S e s and the different approximations to it. The results are given 
in Table [I] for a range of temperatures that includes those relevant to the production of 
8 B in the sun, i.e. 1.3 x 10 7 K < T < 1.6 x 10 7 K |. The values of S cS calculated 
with any of the approximations deviate from the exact results, Eq. (|3]), by no more than 
0.1 eVb or 0.5%. Also, there is less than a 0.1 eVb variation in S e g over the given range 
of temperatures. We note particularly the accuracy of Eq. (|7|), which requires only the 
value of S at 20 keV, and does not require any explicit knowledge of the derivatives or 
of the temperature dependence in E . 

Since the first and second derivatives of the S factor enter into two of the approxima- 
tions, an understanding of the accuracy with which these derivatives can be determined 
is useful. In Table PJ we show the derivatives given in Ref. M for a range of models. 
The first four are hard sphere models [|J depending on the given hard sphere radius, r c . 
The next three are Woods-Saxon models. These models agree with generator coordinate 
model calculations for low energies, E < 300 keV fij, so the conclusions drawn here will 
also apply to results from the generator coordinate model calculations. 

The first derivative shows more model dependence than the second derivative. How- 
ever even for the first derivative the model dependence is quite small. Neglecting the 
hard sphere model with r c = 4.1 fm, which was shown to be unrealistic |3J, and the 
results from Adelberger et al. |3[], the average of the values of the first derivative is 
(1/ S)(dS / cIE)\e=o = — 1.87±0.1 MeV -1 . This introduces a variation of 0.2% in S e s and 
is negligible. 

There is surprisingly little variation in the second derivative. In the hard sphere model 
it does not depend on the radius over the range shown. The results calculated using the 
Woods-Saxon models show only a slight variation. Part of this variation may be due 
to the difficulty in determining the second derivative from the numerical calculations. 
We recommend [l/(2S)](d 2 S/dE 2 )\ E=0 = 15.7 MeV -2 ; the uncertainty in the second 
derivative introduces a negligible error in S e g. 

There has been considerable disagreement in the literature on what values should be 
used for the derivatives. (See for example, Ref. 0). The last line of Table |I| gives values 
of the derivatives from Adelberger et al. |§ that are very different from those we have 
obtained in any of the other models. However, we can reproduce their results by following 
their procedure for determining the derivatives and fitting S(E) by a quadratic form over 
the range E = 20 - 300 keV. Similarly, we can reproduce the values Barker obtained || 
for the derivatives by fitting a quadratic form over the range he chose, - 100 keV. An 
accurate determination of the derivatives from a fit of a quadratic form could only be 
obtained if the fit region was restricted to E < 10 keV. Thus the differences in the values 
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of derivatives quoted in the literature are due to the range of energies used for the fits 
and are unrelated to the true model dependence of the derivatives at threshold. 

We show in Table [TTT] the results of using Eq. (0) with two different sets of derivatives. 
The first is that presented above, namely —1.87 MeV -1 and 15.7 MeV -2 , and the second 
is that from Adelberger et al. given in Table [TT]. Using this second set of derivatives 
introduces an error of 0.3 eVb in S e R. As the reaction rate must be known to better 
than 5% ||, this 1.5% error is 30% of the total error allowed. Remarkably, using the 
inaccurate values of the derivatives at threshold is worse than neglecting the derivatives 
entirely and using S(20). 

In conclusion the use of approximate forms for S e g does not generate significant errors 
if accurate parametrisations of the S factor and its derivatives are used. Eq. (|7|) is of 
special interest since it avoids the use of the derivatives by using S evaluated at 20 keV. 
We therefore recommend that determinations of the 5* factor for 7 Be(p, 7) 8 B at low energy 
quote S(20) as well as S(0). 
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TABLES 



TABLE I. The values of S e g, in eVb, from the different approximations using Eq. (||) for 
S(E). E is in keV and T is in 10 6 K. 



T 


Eq 


Seflf 


^eff-S 


S'cfr-T 


Soff-Snd 


^-820 


12 


15.4 


18.67 


18.68 


18.69 


18.70 


18.58 


13 


16.3 


18.65 


18.66 


18.68 


18.69 


18.59 


14 


17.1 


18.64 


18.65 


18.66 


18.67 


18.59 


15 


17.9 


18.62 


18.63 


18.64 


18.65 


18.60 


16 


18.7 


18.60 


18.61 


18.63 


18.64 


18.60 


17 


19.4 


18.58 


18.59 


18.61 


18.62 


18.61 



TABLE II. The first and second derivatives of S evaluated at threshold for the models of 
Ref. Q . The nomenclature is also that of Ref . Q . 



Model 




1 dS i . 
SlE < MeV "> 




1 d 2 S w 
28 d& (MeV } 


r c = 0.0 fm 




-1.77 




15.7 


r c = 1.0 fm 




-1.82 




15.7 


r c = 2.4 fm 




-1.92 




15.7 


r c = 4.1 fm 




-2.09 




15.7 


Bl 




-1.91 




16.1 


B2 




-1.84 




16.0 


T 




-1.96 




15.7 


Adelberger et al. || 




-0.70 




1.9 


TABLE III. 5 eff _ 


_t, in eVb, 


as obtained using the two different sets of derivatives discussed 


in the text. The temperature is 


in 10 6 K. 






T 


12 


13 14 


15 


16 17 


S eS -T (eVb) 


18.68 


18.66 18.64 


18.62 


18.61 18.59 


S cfr _T (eVb) 


18.97 


18.96 18.96 


18.95 


19.94 18.94 
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